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1.  INTRODUCTION 


Spin  is  usually  imparted  to  artillery  shell  during  launch  through  the  use  of  a  rifled  guri 
tube.  Correspondingly,  most  shell  configurations  include  a  rotating  band,  which  can  be 
described  as  an  aixisymmetric  step-like  protuberance  on  the  projectile  which  niake.s  sliding 
contact  with  the  rifling,  transforming  most  of  the  resulting  torque  into  a  rolling  inotson  of 
the  shell.  Such  rolling  motion  is  normally  required  to  maintain  stable  atmospheric  flight  for 
projectiles  which  do  not  have  lifting  surfaces.  One  aerodynamic  effect  of  the  rotating  band 
is  an  increase  in  the  total  drag  of  the  shell,  thought  to  be  less  than  5%  (.McCoy  1981 ). 

A  protuberance  is  defined  here  to  be  any  physical  irregularity  of  the  projectile  geometry 
(excluding  the  base)  which  has  a  noticeable  influence  on  the  flow.  The  projectile  rotating 
band  can  be  viewed  as  belonging  to  an  extensive  class  of  axisymrnetric  protuberances  which 
are  pertinent  to  many  spinning  and  non-spinning  projectile  configurations.  Other  examples  of 
axisymrnetric  protuberances  include  sabot  grooves,  buttress  threads  or  grcKjve.s,  and  surface 
roughness  in  general. 

The  numerical  prediction  of  the  viscous  flow  over  protuberances  is  an  area  of  interest 
to  the  Army  Research  Laboratory  (ARL).  The  simulation  of  the  flow  field  through  the 
application  of  large-scale  computational  fluid  dynamics  (CFD)  technique-s  prorntses  to  be 
an  accurate  means  for  predicting  the  aerodynamic  effects.  The  simulation  of  flow  over  a 
rotating  band  is  made  difficult  by  the  presence  of  large  flow  gradients  and  separated  flow 
regions. 

In  this  study,  simulation  is  made  of  the  separated,  turbulent  flow  over  a  step- like  protu¬ 
berance  on  a  non-spinning  body  of  revolution  at  Mach  numbers  2.5,  .3.0.  and  3.5  .  A  fimte- 
difference  numerical  technique  is  used  to  solve  the  thin-layer,  compressible,  mass-averaged, 
Navier-Stokes  equations.  The  numerical  computation  is  performed  on  a  non -orthogonal, 
metric-discontinuous,  body-conforming  grid.  The  computed  surface  pressure  distributions 
are  compared  to  existing  supersonic  wind-tunnel  measurements  at  all  three  Mach  numbers. 
Additional  comparison  is  made  with  two  previous  computations  of  this  case  at  Mach  3,0, 
and  tin  evaluation  is  made  of  the  accuracy  and  suitability  of  the  CFD  approach. 

2.  WIND-TUNNEL  MEASUREMENTS 

Experimental  measurements  have  been  made  of  the  surface  pressure  near  a  step- like 
protuberance  on  a  non-spinning  axisymrnetric  body  in  a  supersonic  blow-down  tunnel  located 
at  the  U.S.  Army  Edgewood  Research  Development  and  Engineering  Center,  Edgewood, 
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Maryland  (Dajiberg  and  Palko  1986).  The  wind-tunnel  model.  Figure  1,  is  a  coue-cyliiider 
configuration.  The  diameter  of  the  cylinder,  d,  which  is  taken  as  the  reference  length,  is  2.54 
cm.  The  length  of  the  model  is  approximately  6.8  calibers  (1  cal  1  reference  length),  'i'he 
haJf-angle  of  the  cone  is  13.1®.  Five  pressure  taps  were  installed  in  a  row  along  a  constant 
circumferential  location  on  the  the  cylinder. 

The  step-like  protuberance,  also  in  Figure  1,  is  a  ring  of  inner  diameter  1.0  cal.  outer 
diameter  1.04  cal,  and  width  0.505  cal.  The  ring  wa.s  designed  with  a  set  sere,  in  order  to 
fasten  it  to  a  keyway  inscribed  along  the  cylinder.  The  ring,  positioned  approximately  5  cal 
from  the  nosetip,  could  be  moved  relative  to  the  pressure  taps  so  that  measurements  could 
be  made  at  a  large  number  of  locations  near  the  band.  The  pressure  taps  were  located  on 
the  opposite  side  of  the  cylinder  from  the  keyway.  Measurements  were  made  both  with  and 
without  the  protuberance. 

The  wind-tunnel  model  was  supported  at  the  base  by  an  ajcial  support  and  held  fixed 
relative  to  the  free  stream.  All  measurements  were  made  with  the  model  positioned  at  0° 
incidence.  The  stagnation  temperature  was  260  K.  The  average  Reynolds  number  based 
on  model  diameter  was  1.23,  1.33,  and  1.42  million  at  Mach  2.5,  3.0,  and  3.6,  respectively. 
Data  were  taken  only  after  the  initial  transient  had  decayed,  implying  that,  since  the  model 
was  thin-walled,  the  wall  temperature  was  nearly  adiabatic.  Boundary  layer  transition  was 
initiated  through  the  use  of  a  trip  ring  installed  1  cal  from  the  nosetip,  although  the  entire 
flow  field  was  assumed  to  be  fully  turbulent  in  the  computational  simulations.  A  complete 
tabulation  of  the  experimental  data  is  given  in  the  Appendix. 

Figure  2  is  a  reproduction  of  a  Schlieren  photograph  of  the  cylinder  portion  of  the  wind- 
tunnel  model,  including  the  protuberance,  at  Mach  3.  The  density  gradients  are  exposed  as 
contrasting  shades  of  light.  The  compression  waves  coalesce  into  weak  shock  waves  which 
emanate  from  the  forward  and  rearward  faces  of  the  protuberance.  An  additional  series  of 
compression  waves  upstream  of  the  protuberance  emanate  from  the  keyway  which  secures 
the  ring  to  the  cylinder,  contaminating  a  portion  of  the  flow.  The  pressure  taps  are  located 
on  the  opposite  side  of  the  cylinder,  though,  and  should  not  be  affected. 
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3.  COMPUTATIONAL  APPROACH 


An  unsteady  Navier-Stokes  (UNS)  technique  is  used  to  compute  the  flow  in  the  vicinity 
of  the  rotating  band.  The  UNS  technique  solves  the  governing  equations  by  integrating 
through  time  toward  a  steady-state  solution.  Inflow  conditions  upstream  of  the  rotating 
band  are  provided  by  a  separate  computation  which  utilizes  a  parabolized  Navier-Stokes 
(PNS)  technique.  The  PNS  technique  computes  the  flow  by  integrating  (marching)  in  the 
streamwise  direction.  A  brief  description  of  both  is  included  here,  and  the  references  may 
be  consulted  for  for  additional  details. 

The  use  of  the  PNS  and  UNS  techniques  in  conjunction  is  an  efficient  way  to  n.odel  the 
flow.  The  PNS  technique,  which  processes  only  two  physical  planes  of  data  simultaneously,  is 
applied  over  the  cone  and  cylinder  upstream  of  the  rotating  band,  where  the  flow  is  attached. 
The  UNS  technique,  which  is  more  computationally  intensive,  is  applied  over  the  rotating 
band,  where  separated  flow  regions  exist. 

3.1  Parabolized  Navier-Stokes  (PNS).  Computational  modeling  of  the  flow  over 
the  cone-cylinder  is  accomplished  using  the  PNS  technique  first  reported  by  Schiff  and  Steger 
(1979).  This  technique  spatially  integrates  the  dimensionless,  transformed,  steady,  thin-layer, 
m^lss- averaged  Navier-Stokes  equations  in  strong  conservation  law  form,  given  as: 

-h  /’ij  +  =  Re  ( 1 ) 

Equation  (1)  represents  steady-state  conservation  of  mass,  momentum,  and  energy'  in 
transformed  coordinates  for  large  Reynolds  number  flows.  E,  ,  F,  and  G  are  the  transformed 
inviscid  flux  vectors.  E,  is  a  modified  flux  vector  resulting  from  the  subsonic  sublayer 
approximation.  The  vector  S  is  the  transformed  vector  of  viscous  terms  that  results  from 
the  thin-layer  approximation.  Perfect  gas  behavior  is  assumed.  The  transformed  coordinates 
are:  (  =  ^(x),  the  axial  (marching)  coordinate;  j}  =  T){x,y,z),  the  circumferential  coordinate; 
and  (  =  the  radial  coordinate. 

The  vector  of  dependent  variables  is  defined  as 

Q  =  [p,pu,pv,pw,ef  (2) 

where  the  density  is  p;  the  velocity  components  in  the  x,  y,  and  z  directions  are  u,  v, 
and  lu,  respectively;  and  the  total  energy  per  unit  volume  is  e.  The  solution  is  obtained 
at  each  grid  point  using  the  approximately  factored,  implicit,  delta-form,  finite-difference 
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algorithm  of  Beam  and  Warming  (1978).  Second-order  central-differencing  is  used  in  the 
circumferentiaJ  and  radial  directions,  and  first-order  one-sided  differencing  is  used  in  the 
marching  direction.  Fourth-order  explicit  smoothing  is  added  to  suppress  high  frequency 
oscillations  in  the  solution.  The  solution  is  advanced  downstream  by  numerically  integrating 
in  the  direction  parallel  to  the  projectile  axis.  Each  marching  step  requires  a  series  of  block 
tridiagonal  matrix  inversions  (sweeps)  in  the  circumferential  and  radial  directions. 

At  the  body  surface,  the  no-slip  condition  is  enforced  and  the  pressure  is  obtained  from 
the  subsonic  sublayer  approximation  (i.e.,  the  pressure  across  the  subsonic  portion  of  the 
viscous  layer  is  held  constant).  The  outer  boundary,  which  consists  of  the  bow  shock,  is 
shock  fitted  using  an  implicit  procedure  (Rai  and  Chaussee  1983).  Initial  conditions  for 
marching  are  generated  using  the  PNS  method  in  step-back  mode  (Rai  and  Chaussee  1983; 
Schiff  and  Sturek  1981;  Weinacht  et  al.  1985),  which  assumes  conical  flow  conditions  near 
the  nosetip  and  iteratively  refines  the  solution  to  satisfy  this  assumption. 

The  PNS  computations  provide  inflow  conditions  for  the  ensuing  UNS  computations. 
The  dependent  flow  variables  are  linearly  interpolated  onto  the  UNS  computational  grid 
15  step  heights  upstream  of  the  protuberance.  Additionally,  the  PNS  method  was  used 
in  earlier  portions  of  this  study  to  verify  the  integrity  of  the  measured  pressures  without 
a  protuberance  (Danberg  and  Palko  1986),  and  to  attempt  a  simplified  ramp  model  to 
approximate  the  rotating  band  flowfield  (Danberg  1983). 

Turbulence  is  accounted  for  in  both  the  PNS  and  UNS  computations  using  the  two-layer, 
algebraic  eddy  viscosity  model  of  Baldwin  and  Lomax  (1978).  The  inner  wall-layer  eddy 
viscosity  is  computed  using  a  conventional  Prandit  mixing  length  and  Van  Driest  damping. 
The  outer  or  wake  layer  viscosity  is  based  on  evaluating  the  maximum  in  the  moment  of 
vorticity  and  its  normal  distance  from  the  wall.  The  calculation  of  the  eddy  viscosity  is 
lagged  by  one  time  step  in  the  UNS  code  and  by  one  marching  step  in  the  PNS  code. 

3.2  Unsteady  Navier-Stokes  (UNS).  Computational  modeling  of  the  flow  over  the 
rotating  band  is  accomplished  using  the  UNS  approach  first  reported  by  Pulliam  and  Ste- 
ger  (1978).  This  technique  temporally  integrates  the  dimensionless,  transformed,  time- 
dependent,  thin-layer,  mass- averaged,  Navier-Stokes  equations  in  strong  conservation  law 
form,  given  as: 


Qr  -f  -f  -f  =  Re 


(3) 


Equation  (3)  represents  conservation  of  mass,  momentum,  and  energy  in  transformed 
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coordinates  for  large  Reynolds  number  flow.  The  vector  Q  is  the  transformed  vector  of  de¬ 
pendent  variables.  The  vectors  E,  F,  and  G  are  the  transformed  inviscid  flu.v  vectors.  The 
vector  5  is  the  transformed  vector  of  viscous  terms  that  result  from  the  thiu-layer  approx 
imation  for  large  Reynolds  number  flows.  Ideal  gas  behavior  is  assumed.  1  he  transformed 
coordinates  are  ;  t  =  t,  time;  ^  z),  the  longitudinal  (streamwise)  coordinate;  t)  — 

r}{t,x,y,z)y  the  circumferential  coordinate;  and  (  =  ({t,x,y,2),  the  near-normal  coordinate. 

The  vector  of  dependent  variables,  given  in  Eq.  (2),  is  obtained  at  each  grid  point  using 
the  approximately  factored,  implicit,  delta-form,  finite-difference  algorithm  of  Beam  and 
Warming  (1978).  Second-order  central-differencing  is  used  in  the  three  coordinate  directions, 
and  first-order,  one-sided  time  differencing  is  used.  The  solution  at  each  time  step  requires  a 
series  of  block-tridiagonal  matrix  inversions  (sweeps)  in  each  of  the  transformed  coordinate 
directions.  Fourth-order  explicit  and  second-order  implicit  smoothing  terms  are  added  to 
suppress  high  frequency  oscillations  in  the  solution. 

On  the  body  surface,  the  no-slip  boundary  condition  is  enforced.  The  surface  pressure  is 
found  by  assuming  a  zero  normal  gradient  at  the  wall.  The  density  and  energy  are  found  by 
assuming  adiabatic  wall  conditions.  For  a  constant  coefficient  of  thermal  conductivity,  this 
condition  translates  into  a  zero  density  gradient  in  the  normal  direction,  and  was  applied 
as  a  first-order  approximation  in  the  C  direction.  On  the  outer  boundary,  the  dependent 
flow  variables  are  assumed  free  stream.  On  the  outflow  boundary,  the  dependent  flow  vari¬ 
ables  are  determined  using  linear  extrapolation  along  constant  y  and  grid  lines,  which 
corresponds  to  a  constant  outflow  gradient  boundary  condition.  Axisymmetry  is  invoked  ap¬ 
plying  the  procedure  reported  by  Rizk,  Chaussee,  and  McRae  (1982),  which  employs  three 
circumferential  grid  planes.  Turbulence  is  accounted  for  using  the  algebraic  eddy  viscosity 
model  of  Baldwin  and  Lomax  (1978). 

3.3  Flowfield  Initialization  for  UNS.  Flowfield  initialization  proved  to  be  an  im¬ 
portant  consideration  for  ensuring  numerical  stability  in  the  early  stages  of  the  time-dependent 
computations.  It  is  well-known  that  the  best  initial  flow  field  is  one  which  somewhat  resem¬ 
bles  the  final,  converged  solution.  The  first  initialization  scheme  attempted  here  (beyond  a 
simple  free-stream  initialization,  which  was  inadequate)  was  of  the  form 

Qij,k  =  Qi,j,k  (4) 

where  the  indices  i,j,  and  k  represent  grid  nodes  in  the  and  ^  directions,  respectively. 
That  is,  the  flow  variables  at  each  k  (outward)  index  were  assigned  values  from  the  corre¬ 
sponding  k  index  at  the  inflow  boundary.  This  initialization  was  found  to  be  inadequate 
because  the  computed  solution  quickly  formed  a  large  flow  discontinuity  along  the  forward 
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face  of  the  protuberance,  resulting  in  numerical  instability. 


To  overcome  this  problem,  the  initial  u  velocity  components  of  the  the  first  flowfield 
initialization  were  made  more  incipient,  or  ready  to  separate,  by  multiplying  u  at  each  grid 
point  downstream  of  the  inflow  boundary  by  a  factor  /*,  i.e., 

{'^i,3,k)new  “  (^i,j,k)oidfk  ('^) 


where 


f  ^ 

I  ^mid  1 


for  k  < 


(6) 


[  1.0  for  k  >  kmtd 

The  index  k^id  was  set  equal  to  approximately  half  of  Armor-  This  re-initialization  procedure 
W2is  sufficient  to  allov^  the  solution  to  remain  numerically  stable. 


3.4  Other  Numerical  Aspects.  The  numerical  computations  were  performed  on  a 
Cray  X-MP/48  computer  located  at  the  Ballistic  Research  Laboratory.  Each  computation 
required  322  thousand  words,  or  2.576  million  bytes,  of  computer  memory.  The  computer 
speed  was  approximately  3  seconds  per  integration  step  using  a  140  x  65  computational 
grid.  Approximately  3  hours  of  computer  time  were  used  to  advance  each  numerical  solution 
through  3500  integration  steps.  Additional  integration  steps  were  made  but  did  not  change 
the  solutions  appreciably. 

The  time  increment,  At,  was  determined  by  specifying  a  maximum  allowable  Courant 
number,  Ct,  in  the  computational  domain  (Kutler,  Chaussee,  and  Lombard  1978).  The 
initial  value  of  Ct  was  set  at  5.0  initially,  and  slowly  raised  to  a  value  of  20.0  by  3500  time 
steps.  The  values  of  the  smoothing  parameters  c,  and  (implicit  and  explicit,  respectively) 
were  set  at  0.1  and  0.05  initially,  and  lowered  to  0.06  and  0.03  after  3500  time  steps.  No 
attempt  was  made  to  examine  the  effect  of  various  combinations  of  time  step  and  smoothing 
parameters. 


4.  COMPUTATIONAL  GRIDS 

A  major  point  of  interest  is  the  applicability  of  a  grid  that  “wraps  around”  the  step-like 
protuberance.  For  this  problem,  such  an  approach  requires  the  construction  of  a  highly 
skewed  grid  in  regions  of  sharp  surface  discontinuities. 

The  grid  is  generated  using  a  technique  based  on  the  well-known  work  of  Thompson, 
Thames,  and  Mastin  (1974).  The  governing  equations  are  the  Poisson  equations  in  two 
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dimensions: 

6*  +  ^vv  =  ^(6C) 

Gx  +  Cvv  =  CU.C)  (7) 

In  this  set  of  elliptic  equations,  (  =  i{x,y)  and  (  =  C(^.y)  are  the  transformed  coordinates. 
The  inhomogeneous  terms  P  and  Q  are  forcing  functions  which  control  the  point  spacing 
throughout  the  interior  of  the  domain. 

Steger  and  Sorenson  (1979)  implemented  new  boundary  conditions  to  the  governing  equa¬ 
tions  to  achieve  two  effects:  (1)  to  cluster  points  at  a  specified  distance  from  a  boundary,  and 
(2)  to  force  grid  lines  to  intersect  the  boundary  at  a  specified  angle.  Their  approach  restricts 
the  application  of  the  forcing  terms  P  and  Q  to  the  wall  boundary  ((  =  Cj  ).  which  somewhat 
simplifies  the  equations  and  adds  additional  control  to  the  grid  at  the  wall  boundary.  Their 
condition  that  grid  lines  intersect  the  inner  boundary  in  a  specified  fashion  is  obtained  using 
the  relation 

V^VC=IV^  II  VCl  COS0  (8) 

where  V  is  the  gradient  and  6  is  the  orthogonality  angle,  defined  here  as  the  angle  between  a 
constant  C  Hoe  and  a  constant  f  line.  Using  the  transform  relations  =  y(/J,  etc.,  where 
J  is  the  Jacobian  of  the  transformation),  this  condition  becomes 

+  y«y<  =  -  ({xj  +  yj)(x|  +  y^*)]  cosB  (9) 

where  0  is  specified  along  the  wall  boundary. 

The  grid  generation  scheme  was  modified  here  to  control  0  along  the  wall  by  defining  it 
as  a  function  of  x/d.  To  allow  such  a  function  to  exist  along  the  vertical  faces  of  the  rotating 
band,  the  bottom  corner  of  the  forward  face  of  the  protuberance  was  shifted  upstream  by 
0.001  cal,  from  x/d=4.95  to  x/d— 4.949.  The  change  in  vertical  orientation  of  the  forward 
face  due  to  this  shift  is  approximately  1.43'*,  and  is  not  expected  to  change  the  flow  field 
from  that  of  a  vertical  face  (Young  and  Paterson  1981).  The  orthogonality  angle  along  the 
wall  was  prescribed  upstream  of  the  midpoint  of  the  band;  downstream  of  the  midpoint,  the 
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grid  was  constructed  using  reflection  principles.  The  orthogonality  angle  was  prescribed  as: 

for  4.35  <  xfd  <  4.949  (upstream  of  band) 

for  xfd  —  4.949  (bottom  front  corner  grid  point) 

for  4.949  <  xjd  <  4.95  (front  face  of  band)  (10) 

for  xjd  =  4.95  (top  front  corner  grid  point) 

for  4.95  <  xfd  <  5.20  (top  front  half  of  band) 

The  computational  grid  used  for  the  Mach  3.0  c<ise  is  shown  in  Figure  3.  In  the  longitudi¬ 
nal  direction,  140  grid  points  were  used.  In  the  outward  direction,  65  grid  points  were  used. 
The  grid  does  not  change  with  time.  The  inflow  boundary  is  located  at  r/d=4.35,  or  15  step 
heights  upstream  of  the  forward  face  of  the  protuberance.  The  outflow  boundary  is  located 
at  a:/d=6.05,  or  15  step  heights  downstream  of  the  rearward  face  of  the  protuberance. 

Figures  4  and  5  show  closeup  views  of  the  computational  grid  at  the  forward  and  rearward 
faces  of  the  protuberance.  As  usual  for  viscous  computations,  the  grid  is  clustered  near  the 
wall.  The  distance  between  the  wall  and  the  first  outward  grid  point  is  held  constant  for 
each  case.  The  values  were  obtained  from  the  corresponding  PNS  solution  at  the  inflow 
plane,  which,  through  an  adaptive  grid  technique  (Sturek  and  Schiff  1981),  maintained  the 
outward  spacing  at  the  wall  such  that  the  first  grid  point  from  the  wall  satisfied  the  condition 
3  <  <  5.  The  longitudinal  point  distribution  is  prescribed  such  that  grid  points  were 

slightly  clustered  along  the  forward  and  rearward  faces.  Each  face  contains  9  grid  points, 
including  the  top  and  bottom  corners. 

Appreciable  discontinuities  in  the  grid  metrics  extend  outward  from  the  top  and  bot¬ 
tom  corners  of  each  face.  The  grid  is  highly  skewed  in  these  regions,  at  about  45°.  Grid 
non-orthogonality  at  the  wall  is  an  issue  of  interest  since  the  wall  boundary  conditions  are 
implemented  with  the  assumption  that  the  grid  is  orthogonal  to  the  wall.  The  effective  error 
is  expected  to  be  small  due  to  small  grid  spacing  in  the  outward  direction  as  well  as  the 
small  longitudinal  gradients  (compared  to  the  outward  gradients)  which  occur  at  the  wall. 


B={ 


2  ~  4  (4.494  -  4.35) 

JT 
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3t 
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5.  RESULTS 


5.1  Boundary  Layer  Dispiacement  Thickness.  The  relative  sizes  of  the  oucom- 
ing  boundary  layers  with  respect  to  the  rotating  band  thickness  give  an  indication  of  the  scale 
of  the  problem  (Sedney  1973).  The  boundary  layer  displacement  thickness,  (5*,  was  deter¬ 
mined  from  the  converged  numerical  solutions  15  step-hights  upstream  of  the  rotating  band. 
A  compressible,  axisymmetric  boundary  layer  was  assumed;  the  displacement  thickness  was 
calculated  using 

8^ 


fiy- 


pu 


Ptut 


)rdr  +  r^-r, 


in) 


where  the  subscripts  w  and  c  represent  wall  and  edge  conditions,  respectively.  At  each 
Mach  number,  a  variation  of  8”  was  generated  as  a  function  of  assumed  edge  location.  The 
Ue  profiles  were  inspected  and  found  to  be  well-behaved.  The  displacement  thickness  was 
maximum  near  the  apparent  boundary  layer  edge  and  the  edge  location  was  assumed  to 
occur  at  this  maximum.  The  accurjicy  of  the  predicted  6*  is  within  1%;  the  determined 
values  are  shown  in  Table  1.  The  table  shows  that,  in  each  case,  the  computed  displacement 
thickness  is  about  half  of  the  rotating  bcind  thickness. 


Table  1.  Computed  Boundary  Layer  Dis  acement  Thicknesses  at  x/d=4.35  . 


M 

8"  (cal) 

2.5 

0.0159  ±  1% 

3.0 

0.0186  ±  1% 

3.5 

0.0226  ±  1% 

5.2  Flow  Structure.  Figures  6  and  7  show  computed  velocity  vectors  for  the  Mach  3 
case  near  the  forward  and  rearward  faces  of  the  rotating  band,  respectively.  Longitudinally- 
separated  flow  regions  are  captured  in  the  numerical  solution  both  upstream  and  down¬ 
stream  of  the  band.  The  computed  separation  and  reattachment  points  were  more  precisely 
determined  at  each  Mach  number  by  examining  the  computed  skin  friction  coefficient.  The 
velocity  derivative  used  to  compute  the  skin  friction  was  approximated  using  a  one-sided, 
first-order  difference.  The  skin  friction  was  assumed  to  vary  linearly  between  adjacent  grid 
points;  separation  or  reattachment  was  assumed  to  occur  where  the  skin  friction  was  zero. 
Table  2  shows  the  computed  separation  and  reattachment  distances  for  each  of  the  three  dif¬ 
ferent  Mach  numbers.  The  separation  distance  is  computed  to  vary  by  less  than  1%  between 
the  three  Mach  numbers.  The  non-monotonic  variation  of  separation  distance  with  respect 
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to  Maoh  number  is  likely  within  the  accuracy  of  the  computation,  given  the  current  longi¬ 
tudinal  grid  resolution  and  the  assumptions  used  to  compute  skin  friction  and  separation 
location. 


Table  2.  Computed  Separation  and  Reattachment  Distances  (in  Step  Heiglits). 


M 

(Sa  Sfgp)/h 

(■Sreat  ^b)/^ 

2.5 

5.05 

3.99 

3.0 

5.03 

3.59 

3.5 

5.09 

3.37 

Figures  8  2ind  9  show  the  computed  Mach  contours  in  the  forward  and  backward  step 
regions  for  the  Mach  3.0  case.  The  sonic  line  is  emphasized.  Only  slight  jaggedness  in  the 
Mach  contours  are  discernable  along  the  constant  ^  grid  lines  which  emanate  from  the  top 
and  bottom  corners  of  the  faces,  where  large  metric  discontinuities  exist.  Figure  10  shows 
computed  Mach  contours  for  the  entire  computational  domain.  The  crispness  of  the  weak 
shocks  in  the  far  field  is  definitely  limited  by  the  poor  grid  resolution  there.  The  compression 
which  emanates  from  the  front  face  is  consistently  captured  between  three  adjacent  grid 
points,  but  far  from  the  body  such  capturing  still  gives  a  largely  smeared  structure.  The 
recompression  just  downstream  of  the  rearward  face  appears  to  be  captured  slightly  better, 
possibly  because  of  better  grid  resolution  there.  The  computed  bow  shock  is  also  smeared 
downstream  from  the  inflow  boundary,  and  the  disturbances  propagate  toward  the  body  and 
interact  with  the  waves  which  emanate  from  the  rotating  band.  The  influence  of  these  outer 
disturbances  on  the  computed  flow  field  near  the  rotating  band  is  thought  to  be  negligible: 
first,  the  characteristic  lines  are  not  of  sufficient  slope  to  intersect  the  body  close  to  the  band, 
and  second,  the  flow  is  nearly  uniform  in  the  outer  region. 

Figures  11  and  12  show  the  computed  pressure  contours  in  the  forward  and  backward 
step  regions  for  the  Mach  3.0  case.  The  contours  are  plotted  as  p/Poo-  The  highest  surface 
pressure  on  the  front  face,  occurring  near  the  top  of  the  band,  is  computed  to  be  2.91. 
The  stagnation  region  appears  to  be  spread  over  several  grid  points,  rather  than  being 
concentrated  at  a  particular  point.  Some  jaggedness  is  present  in  the  contours  near  the 
top  corner  of  the  forward  face.  The  lowest  surface  pressure  on  the  rearward  face,  occurring 
near  the  top  of  the  band,  is  computed  to  be  0.104.  The  computed  expansion  begins  at 
about  two  grid  points  upstream  of  the  corner  of  the  rearward  face  (x/d  «  5.435),  indicating 
that  flow  disturbances  are  being  propagated  upstream  through  the  subsonic  portion  of  the 
boundary  layer.  Jaggedness  in  the  contours  near  the  top  corner  of  the  rearward  face  is  not 
nearly  as  apparent  as  that  near  the  forward  face,  possibly  because  the  gradients  are  not  as 
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laxge.  Figure  13  shows  computed  pressure  contours  for  the  entire  computational  domain. 
The  contour  which  runs  nearly  parallel  and  next  to  the  outer  boundary  is  the  level  l  .i,  and 
represents  the  captured  outer  bow  shock. 


5.3  Surface  Pressure  Comparison.  The  pressure  is  compared  in  the  form  ACp, 
which  is  the  difference  in  surface  pressure  coefficient  with  and  without  the  rotating  band, 


i.e., 


Aa 


'PI 


Pl/Poo  -  Po/Po 


(12) 


where  pi  and  po  are  the  surface  pressures  with  and  without  the  band,  respectively.  The 
formal  accuracy  of  the  pressure  measurements  is  discussed  in  detail  in  the  experiment  doc¬ 
umentation  (Danberg  and  Palko  1986).  The  repeatability  of  the  measured  ACp  appears  to 
be  approximately  0.01  . 


The  surface  pressure  distributions  are  plotted  as  functions  of  distance  along  the  surface 
coordinate,  s  (which  is  coincident  with  ^),  scaled  by  the  band  height,  h.  The  origin  of  s  for 
comparisons  upstream  of  the  band  is  Sa,  located  at  the  bottom  corner  of  the  forward  face. 
The  origin  of  a  for  comparisons  downstream  of  the  band  is  s^,  located  at  the  bottom  corner 
of  the  rearward  face.  The  final,  tabulated  data  is  listed  in  the  Appendix. 

Upstream  of  the  rotating  band.  Figures  14-16,  the  computed  and  measured  pressures 
compare  to  within  about  0.01  in  the  attached  flow  region  at  all  three  Mach  numbers.  In  all 
three  cases,  the  location  of  initial  pressure  rise  associated  with  the  onset  of  flow  separation 
is  predicted  to  occur  further  upstream  than  the  measurements  show.  Near  this  location,  the 
computed  ACp  is  as  much  as  0.03  greater  than  the  measured  value.  At  about  2  to  4  step 
heights  upstream  of  the  band,  the  difference  in  the  computed  and  measured  ACp  is  within 
0.01.  At  less  than  2  step  heights  upstream  of  the  band,  the  computation  overpredicts  ACp  to 
a  maximum  of  about  0.04  at  Mach  2.5  and  3.0.  At  Mach  3.5,  the  computation  overpredicts 
ACp  to  a  maximum  of  about  0.02. 

Downstream  of  the  rotating  band,  Figures  17-19,  the  computed  and  measured  pressures 
compare  quite  differently  than  upstream  of  the  band.  Within  1  step  height  downstream  of 
the  band,  the  computed  and  measured  pressures  agree  within  about  0.015.  However,  the 
experiment  shows  a  rapid  pressure  recovery  2  to  4  step  heights  downstream  of  the  face, 
whereas  the  computation  shows  a  slow  pressure  recovery  extending  beyond  8  step  heights 
downstream.  This  trend  is  repeated  at  all  three  Mach  numbers. 


Figures  15  and  18  include  the  results  of  a  previous  computation  of  the  Mach  3.0  case  (Sahu 
1986).  Those  results  were  generated  using  the  same  numerical  scheme  as  in  this  study,  but 
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with  some  important  differences  in  the  approach.  The  first  difference  is  that  the  computation 
was  carried  out  on  an  orthogonal  grid,  as  opposed  to  the  highly  skewed  grid  of  this  study. 
Second,  the  no-sIip  boundary  condition  was  applied  along  the  forward  and  rearward  step 
faces  by  using  a  fiowfield  blanking  procedure.  Third,  the  Baldwin-Lomax  turbulence  model 
was  only  applied  along  grid  lines  which  emanated  from  the  cylinder  surface,  meaning  that 
the  effective  turbulence  was  not  adequately  specified  along  the  forward  and  rearward  step 
faces.  Fourth,  the  longitudinal  grid  resolution  in  the  separated  flow  regions  was  not  as  good 
as  in  this  study.  In  Figures  15  and  1 8,  the  two  computational  results  are  similar  in  that 
they  agree  with  some,  but  not  all,  of  the  measured  pressures  upstream  of  the  band.  Close 
to  the  forward  face,  the  prediction  of  the  present  study  appears  to  agree  with  the  measured 
pressure  slightly  better  than  that  of  Sahu.  Atop  of  the  band,  the  computations  disagree  with 
each  other,  but  the  prediction  of  the  present  study  appears  to  agree  better  with  the  single 
available  data  point.  Downstream  of  the  band,  both  computations  fail  to  predict  the  rapid 
pressure  recovery  that  occurs  from  1  to  3  step-heights  downstream.  Sahu’s  computation 
shows  a  faster  pressure  recovery  than  the  present  computation,  but  this  is  likely  due  to  its 
higher  pressure  atop  of  the  band  itself.  No  comparison  can  be  made  on  the  vertical  faces  of 
the  band,  since  Sahu’s  values  are  no  longer  available  and  were  never  published. 

Still  another  computation  of  the  Mach  3.0  case  has  been  reported  (Danberg  and  Patel 
1989).  Those  results  were  generated  by  solving  the  full,  mass-averaged  Navier-Stokes  equa¬ 
tions  using  a  zonal,  explicit  numerical  approach  on  an  orthogonal  grid.  The  Baldwin-Lomax 
turbulence  model  was  used  in  a  modified  form.  Through  an  examination  of  the  reported 
pressure  distribution,  it  is  easily  determined  that  agreement  with  the  measured  surface  pres¬ 
sure  was  quite  similar  to  the  present  computation  and  to  that  of  Sahu.  Upstream  of  the 
band,  the  pressure  near  the  forward  face  is  overpredicted.  Downstream  of  the  band,  the 
rapid  pressure  recovery  is  not  predicted. 

The  three  independent  computational  results  suggest  that  the  thin-layer  model  and  the 
highly  skewed  grid  are  not  major  sources  of  error  in  the  present  study.  The  experimental 
measurements  are  assumed  to  be  accurate.  It  is  concluded  the  turbulence  model  is  the  most 
likely  contributor  to  inaccur<icies  in  the  computed  pressure  field  in  this  study,  as  well  as  in 
the  other  two  studies. 

The  inaccuracy  of  the  computed  recovering  pressure  field  downstream  of  the  band  raises 
some  concern  for  the  applicability  of  the  CFD  approach.  The  current  CFD  capability  may 
have  value  for  predicting  the  drag  due  to  the  presence  of  the  rotating  band,  although  confir¬ 
mation  of  this  has  not  yet  been  made.  Experimental  drag  data  with  and  without  a  rotating 
band  would  be  valuable.  The  computed  recovering  boundary  layer  is  probably  no  more 
accurate  than  the  computed  recovering  pressure  field.  Computations  of,  for  instance,  Mag- 


nus  effects  (which  are  critically  dependent  on  the  boundary  layer  profile)  would  need  to  be 
validated.  The  velocity  data  which  has  been  published  by  Gorney,  Yanta,  and  Ausherman 
(1982)  may  be  valuable  for  validating  the  accuracy  of  the  computed  attached,  recovering 
boundary  layer  downstream  of  the  band. 


6.  CONCLUSION 

The  numerical  solution  showed  the  displacement  thickness  of  the  oncoming  boundary- 
layer  to  be  the  same  order  of  magnitude  as  the  step  height.  Regions  of  flow  separation  were 
captured  in  the  numerical  solutions.  In  the  recirculation  region  upstream  of  the  protuber- 
tince,  agreement  in  computed  and  measured  pressures  was  within  the  measurement  accuracy 
at  most  tap  locations.  Near  the  predicted  separation  point  and  within  2  step  heights  of  the 
forward  face,  agreement  was  not  within  the  measurement  accuracy.  In  the  interaction  re¬ 
gion  downstream  of  the  protuberance,  agreement  was  within  the  merisurement  accuracy  only 
within  1  step  height  from  the  rearward  face,  and  the  rate  of  pressure  recovery  w<is  underpre¬ 
dicted  in  the  numerical  solutions.  The  pressure  comparison  upstream  of  the  protuberance 
was  significantly  better  than  the  comparison  downstream  of  the  protuberance. 

A  highly  skewed,  metric-discontinuous,  wrap-around  grid  was  used  to  perform  the  compu¬ 
tations.  Except  for  some  minor  jaggedness  in  the  pressure  contours,  no  major  discontinuities 
can  be  discerned  in  the  numerical  solution  where  the  grid  metrics  change  rapidly.  Two 
other  CFD  studies  which  did  not  use  highly  skewed  grids  were  found  to  have  similar  dis¬ 
agreement  in  surface  pressure  as  this  study.  Under  these  considerations,  the  non-orthogonal, 
metric-discontinuous  grid  does  not  appear  to  be  an  obvious  cause  of  soluti  jn  inaccuracies. 
A  more  likely  contributor  to  inaccuracies  in  the  computed  solutions  of  all  three  studies 
is  the  Baldwin-Lomax  algebraic  turbulence  model,  which  has  been  employed  here  without 
modification. 

Further  study  in  the  modeling  of  separated  flow  over  protuberances  will  need  to  address 
the  issue  of  turbulence  modeling  in  detail.  Higher  order  models  may  be  needed  to  improve 
the  accuracy  of  the  numerical  prediction.  Attention  should  be  given  to  the  accuracy  of  the 
reattached,  recovering  boundary  layer  downstream  of  the  rotating  band.  An  experimental 
consideration  of  the  drag  increase  due  to  the  rotating  band  would  be  of  tremendous  value  to 
evaluate  CFD  models.  An  extension  of  the  present  work  could  include  a  substantial  increase 
in  grid  resolution;  the  present  grid  resolution  was  chosen  to  roughly  match  that  of  the  two 
previous  computational  studies  which  were  cited. 
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Ail  Dinensions  in  Calibers 
Dianeter,  d=254cn 


Figure  1.  Iliustration  of  Wind  Tunnel  Model  and  Rotating  Band 
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Figure  2.  Schlieren  Photograph  of  Flow  over  Step-like  Protuberance  at  Mach  i 
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Figure  3.  UNS  Computational  Grid,  A/=3.0 
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Figure  6.  Computed  Velocity  Vector 
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=3.0  (0.1  to  3.2  in  increments  of  0.1) 


Figure  13.  Computed  Pressure  Contours,  Af=3.0  (0.1  to  3.0  in  increments  of  0.1) 
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(s-Sa)/h 


Figure  15.  Surface  Pressure  Comparison  near  Forward  Face,  M=3.0 
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Figure  16.  Surface  Pressure  Comparison  near  Forward  Face,  Af=3.5 
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Figure  17.  Surface  Pressure  Comparison  near  Rearward  Face,  M=2.5 
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Figure  18.  Surface  Pressure  Comparison  near  Rearward  Face,  M=3.0 
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Figure  19.  Surface  Pressure  Comparison  near  Rearward  Face,  Af=3.5 
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APPENDIX: 

TABULATED  SURFACE  PRESSURE  DATA 


_ M==2.5 _ 

upstream  of  band  ||  downstream  of  band 


-14.9 

hubssohii 

0.312 

-0.142 

-15.6 

-0.001 

0.500 

-0.142 

-14.4 

-0.001 

1.125 

-0.144 

-13.9 

-0.004 

1.625 

-0.129 

-12.9 

-0.005 

2.125 

-0.092 

-12.8 

-0.003 

2.375 

-0.075 

-12.4 

-0.009 

2.625 

-0.061 

-12.4 

-0.005 

3.125 

-0.037 

-11.4 

-0.003 

3.625 

-0.029 

-10.5 

-0.002 

4.125 

-0.028 

-9.59 

-0.004 

4.625 

-0.018 

-9.40 

■'H'l'M 

5.125 

-0.023 

-8.77 

5.150 

-0.023 

-8.60 

-0.002 

6.150 

-0.022 

-8.28 

-0.003 

7.025 

-0.018 

-7.97 

mmm 

7.837 

-7.79 

8.025 

-0.014 

-7.78 

-0.002 

8.650 

-0.014 

-7.36 

-0.002 

8.900 

-0.010 

-7.26 

-0.001 

9.150 

-0.012 

-7.17 

-0.002 

9.525 

-0.011 

-6.78 

-0.006 

9.650 

-0.009 

-6.26 

-0.005 

9.713 

-0.010 

-6.04 

-0.005 

10.15 

-0.007 

-5.30 

0.003 

10.65 

-0.006 

-4.78 

11.15 

-0.007 

-4.75 

0.019 

12.65 

-0.006 

-4.50 

0.079 

12.67 

-0.002 

0.138 

13.00 

-3.55 

0.153 

13.67 

-0.001 

0.156 

14.00 

-0.004 

-3.05 

0.177 

14.55 

0.0003 

-2.86 

0.184 

15.50 

-0.004 

-2.00 

0.206 

16.42 

0.0002 

-1.00 

0.219 

17.05 

0.0009 

-0.38 

0.235 

17.12 

-0.006 

-0.19 

0.249 

17.24 

-0.00003 

upstream  of  band 


=3.0 


downstream  of  band 


(s  -  st,)/h 


-0.002 


-0.004 


-0.001 


-0.003 


-0.003 


03 


-0.004 


2 


-0.002 


2 


-0.006 


-0.005 


3 


-0.002 


6 


-0.002 


-0.006 


-0.006 


-0.007 


5 


3 


-0.005 


6 


-0.003 


0.001 


0.005 


0.046 


0.099 


0.118 


0.118 


0.138 


0.138 


0.167 


0.169 


0.177 


0.186 


0.194 


0.313 


-0.095 


-0.097 


-0.095 


1.625 

-0.090 

2.125 

-0.067 

2.375 

-0.062 

2.625 

-0.044 

3.125 

-0.026 

3.625 

-0.016 

4.125 

-0.016 

4.625 

-0.015 

5.125 

-0.013 

5.150 

-0.013 

6.150 

-0.012 

7.025 

-0.012 

7.837 

8.025 

-0.012 

8.650 

-0.011 

8.898 

-0.013 

8.900 

-0.011 

9.150 

9.525 

-0.010 

9.650 

-0.010 

mmsm 

-0.014 

10.15 

10.65 

-0.008 

11.15 

-0.009 

12.65 

-0.005 

12.67 

-0.004 

13.00 

-0.004 

13.67 

-0.0001 

14.00 

-0.001 

14.55 

-0.0004 

15.50 

-0.002 

16.42 

-0.001 

16.42 

mmm 

0.002 

17.13 

M=3.5 


upstream  of  baud 

downstream  of  band 

(s  -  Sa)//» 

^goii 

{s  -  Sb)/h 

-14.9 

-0.002 

0.312 

-0.065 

-14.6 

-0.002 

0.500 

-0.066 

-14.4 

1.125 

-0.067 

-13.9 

1.625 

-0.063 

-12.9 

-0.002 

2.125 

-0.052 

-12.8 

-0.009 

2.375 

-0.046 

-12.4 

-0.003 

2.625 

-0.040 

-12.4 

-0.003 

3.125 

-0.027 

-11.4 

-0.002 

3.625 

-0.014 

-10.5 

-0.001 

4.125 

-0.013 

-9.59 

-0.003 

4.625 

-0.012 

-9.40 

5.125 

-8.78 

-0.002 

5.150 

-0.009 

-0.002 

6.150 

-0.010 

-8.59 

-0.002 

7.025 

-0.011 

7.837 

-7.97 

-0.004 

8.025 

-0.009 

-7.79 

-0.002 

8.650 

-0.009 

-7.78 

-0.001 

8.898 

-0.011 

-7.36 

-0.003 

8.900 

-0.011 

-7.26 

9.150 

-0.008 

-7.17 

-0.002 

9.525 

-6.78 

-0.002 

-0.009 

-0.002 

-0.010 

10.15 

-0.009 

-5.30 

-0.0007 

10.65 

-0.009 

-4.78 

0.003 

11.14 

-0.009 

-4.75 

0.010 

12.65 

mm'ssm 

-4.50 

0.006 

12.67 

-0.007 

-3.75 

0.027 

13.00 

-0.008 

-3.55 

0.066 

13.67 

-0.005 

-3.50 

0.089 

14.00 

-0.006 

-3.05 

0.107 

14.55 

-0.007 

-2.88 

0.106 

15.50 

-2.00 

0.137 

16.42 

-0.003 

-1.01 

0.146 

16.42 

-0.003 

-1.00 

0.142 

17.05 

-0.003 

0.153 

17.13 

-0.003 

-0.18 

0.167 

17.24 

■BZBI 

39 


Intentionally  left  blank. 


40 


LIST  OF  SYMBOLS 


a  speed  of  sound 

Cp  surface  pressure  coefficient 

Ct  Courant  number 

Cp  constant  pressure  specific  heat 

d  reference  diameter  of  body  of  revolution  (equal  to  1  caliber) 

jP,  Ginviscid  flux  vectors  of  transformed  gas  dynamic  equations 
h  rotating  band  thickness  (protuberance  or  step  height) 

i,j,  k  grid  point  indices  associated  with  t/,  C,  respectively 
J  Jacobian  of  coordinate  transformation 

L  reference  length  (equal  to  1  caliber) 

M  free- stream  Mach  number 

p  pressure 

Q  vector  of  dependent  variables  of  gas  dynamic  equations 

Q  vector  of  dependent  variables  of  transformed  gas  dynamic  equations 

Re  Reynolds  number,  /9<»UooL//Xoo 

Re  sonic  Reynolds  number,  PooOooLjpoo 

r  radial  coordinate 

S  viscosity  vector  of  transformed  gas  dynamic  equations 

s  arc  length  in  ^  direction 

t  time 

u,  V,  w  velocity  components  in  ar,  y,  z  directions 
u  velocity  component  close  to  and  parallel  to  wall 

X,  y,  z  physical  Cartesian  coordinates 

y  normal  distance  from  wall 


Greek  Symbols 

7  ratio  of  specific  heats 

6*  boundary  layer  displacement  thickness 

Ci,  Ce  implicit  and  explicit  numerical  damping  parameters 

e  total  energy  per  unit  volume  of  fluid 

p  coefficient  of  molecular  viscosity 

p  density 

r  transformed  coordinate  of  time 

6  angle  of  grid  orthogonalit}' 

rj,  C  transformed  coordinates 


Subscripts 

a  evaluated  at  bottom  corner  of  forward  facing  step 

h  evaluated  at  bottom  corner  of  rearward  facing  step 

c  boundary  layer  edge  condition 

Teat  evaluated  at  flow  reattachment  location 
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sejt  evaluated  at  flow  separation  location 

w  wall  condition 

oo  free-stream  condition 

1/  evaluated  at  / 
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Mechanics 

ATTN:  D.  S.  Dolling 
Austin,  Texas  78712-1055 

1  University  of  Delaware 
Department  of  Mechanical 
Engineering 
ATTN:  L.  Schwartz 
Newark,  DE  19716 


1  MDA  Engineering,  Inc. 
ATTN;  J.  P.  Steinbrenner 
500  E.  Border  St. 

Suit  401 

Arlington,  TX  76010 


University  of  Cincinnati 

Department  of  Aerospace  Engineering 

ATTN:  S.  Rubin 

Mail  Location  70 

Cincinnati,  OH  45221 


1  Los  Alamos  National  Laboratory 
ATTN:  Bill  Hogan 
MS  G770 

Los  Alamos,  NM  87545 

1  Southwest  Research  Institute 
ATTN:  David  L.  Littlefield 
San  Antonio,  TX  78228-0510 

2  University  of  Cedifornia,  Davis 
Department  of  Mechanical 

Engineering 
ATTN:  H.A.  Dwyer 
B.  Meakin 
Davis,  CA  95616 

1  University  of  Maryland 

Department  of  Aerospace  Engr. 
ATTN:  J.  D.  Anderson,  Jr. 
College  Park,  MD  20742 


1  University  of  Florida 
Department  of  Engineering  Sciences 
College  of  Engineering 

ATTN:  C.  C.  Hsu 
Gainesville,  FL  32611 

2  Alliant  Techsystems,  Inc. 

ATTN:  M.  Swenson 

R.  Burretta 

Mail  Station  MN48-3700 
7225  Northland  Dr. 

Brooklyn  Park,  MN  55428 

2  Director 

US  Army  ERDEC 
ATTN:  SCBRD-RTB 
D.  Weber 
F.  Wrede 

Aberdeen  Proving  Ground,  MD  21010-5423 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 


1.  ARL  Report  Number  ARL-TR-102 


March  1993 
Date  of  Repon _ 


Your 


2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.)  _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  Slate,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Inconect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State.  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


Depafttment  of  the  Army 


OFROAL  BUSINESS 


BUSINESS  REPLY  iVL\IL 

RRST  ClASS  Mo  oaoi,  aps,  md 

Postage  <»)li  be  paid  by  addiessee 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 
Aberdeen  Proving  Ground,  MD  21005-5066 


